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Abstract 

Recent research on the dynamics of certain fluid dynamical instabilities shows that 
when there is a slow invariant manifold subject to fast timescale instability the 
dynamics are extremely sensitive to noise. The behaviour of such systems can be 
described in terms of a one-dimensional map, and previous work has shown how 
the effect of noise can be modelled by a simple adjustment to the map. Here we 
undertake an in depth investigation of a particular set of equations, using the meth- 
ods of stochastic integration. We confirm the prediction of the earlier studies that 
the noise becomes important when fj\ lne| = 0(1), where \x is the small timescale 
ratio and e is the noise level. In addition, we present detailed information about the 
statistics of the solution when the noise is a dominant effect; the analytical results 
show excellent agreement with numerical simulations. 
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1. Introduction 

In many circumstances, a low order system of ordinary differential equations (ODEs) 
serves as a useful model for a physical system. A difficulty arises, however, if the 
solutions are noticeably affected by small external noise. This is the case in several 
systems of physical interest sharing the characteristic that their solutions consist of 
alternating slow and fast phases [1,2,3]. In this paper we take as our example the 
following third order system of ODEs describing the resonant interaction of three 
wave modes when one mode is unstable and the other two damped [4,5,6]: 

x = fix — y 2 + 2z 2 — 8z, 

y = y(x-l), (1) 
z = nz + Sx — 2xz. 

When [i is small the character of solutions of (1) is dramatically changed by 
tiny amounts of additive noise - the bifurcation structure with a full gamut of 
periodic orbits and chaotic regions is replaced by a noisily periodic orbit across 
a wide range of parameter values (Figure 1). The quantity /j is the ratio of the 
instability of the unstable mode (its exponential rate of growth in the absence of 
interaction) to the damping rates of the other two modes (assumed equal). [7] 

It is possible to describe the dynamics of (1) in terms of a one-dimensional 
map. As shown in [1], analytical expressions can be obtained for this map by 
assuming that the solutions consist of alternating slow and fast phases, and solving 
approximations to (1) in each phase (Figure 2). In the slow phase the system is 
close to the invariant plane y = and moves slowly from the region where the 
plane is attracting (x < 1) to the region where it is repelling (x > 1). This phase 
is occasionally interrupted by a fast phase which reinjects the system close to the 
attracting part of y = 0. It is in the slow phase that the sensitivity to noise arises. 

Here we use a stochastic differential equation to calculate the adjustment to 
the map of Hughes and Proctor [1] necessary to describe the dynamics of the slow 
phase in the presence of additive white noise. The smallness of ft is responsible for 
the division of the dynamics into two phases and the sensitivity to noise. In this 
paper, all calculations are done to lowest order in [i. 

The map we use is a map of successive turning points of x. In the presence 
of noise the turning point of x which defines the end of the slow phase is a random 
variable, x max . In the heuristic model of Hughes and Proctor [1], the noise deter- 
mines x m ax if lA hie| < 0(1), where e is the r.m.s. noise level, and x max is estimated 
by assuming that \y\ = O(e) at x = 1. The value of x max calculated in this manner 
is O(l) less than the corresponding deterministic value. 

In this paper, we extend the treatment of Hughes and Proctor using a stochas- 
tic differential equation to describe the slow phase. We derive explicit expressions for 
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Figure 1. Bifurcation diagrams with and without noise for u, = 0.01. The graphs 
are obtained from numerical simulation of (1 ) and each dot represents a turning 
point of x. The top graph is obtained with no added noise; the bottom graph 
with very small noise (r.m.s magnitude e = 10~ 10 ) added to the variable y. The 
noiseless bifurcation structure including chaotic regions is replaced by a noisily 
periodic orbit for S < 1.5. 
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Figure 2. The one- dimensional map. (If x = at x = xi then the next value of 
x for which x = is xi .) In the top figure, the solid line is the formula of Hughes 
and Proctor [1] for fi = 0.01 and 5 = 0.5. The dots are obtained from numerical 
solution of (1). The same map is shown below with a Hat top due to noise. The 
two lines in the Bat top are (x m ax) + &x max and (x m ax) -o- Xmax , calculated 
as described in the text with e = 10~ 10 , and the dots are numerical results with 
white noise of magnitude 10~ 10 added to the variable y. 

the probability distribution of x max , and for the condition on lne| which marks 
the transition to the noise-controlled regime. In this regime, the most probable 
value of x max is a function of n\ lne| and the standard deviation of its probability 
distribution is proportional to fx. For small fx, our calculations predict accurately 
the results of numerical solution of (1) with low- level noise added to the variable 
y. The numerical results presented here were obtained using a simple extension 
of the Heun (second order Runge-Kutta) method for integrating ODEs to include 
additive white noise [8, 9]. The increment to y at each step includes a Gaussian 
random variable proportional to e and to the square root of the timestep. 

Slow-fast dynamics similar to those of the three-wave resonance system are 
relevant in other contexts. Our results are presented in a such a way that they can 
be easily generalised. In the appendix, we show how a generic slow phase is reduced 
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to a problem taken from dynamic bifurcation theory, and summarise some results 
for this case. 

2. The slow phase with and without white noise 

The slow phase of (1) is defined as beginning when the following are true: 

y 2 <^ fi, z=- + 0(/j,), and 0(n)<x<l. (2) 
We then find that (1) reduce to: 

x = nf(x) -y 2 , 



y = yg(x), (3) 



5 s 



2 



where f(x) = x + — and g(x) = x — 1. (4) 

When the initial conditions (xo, yo, zq) satisfy (2) we observe the following: 
-the variable y decreases exponentially until x = 1 and then increases exponentially; 
-the variable x is the driving variable, evolving independently until the very end of 
the slow phase; 

-the remaining variable, z, is of secondary importance for x > 0(/x) because it is 
'slaved' to x (given function of x). 

Our strategy for determining x max is to take x = fif(x), so that x is a function 
of time, and then solve for y as a function of time. x max is then the value of x at 
which y 2 = fif(x). The term — y 2 in the equation for x has a small effect at the end 
of the slow phase which we calculate a correction for. 

In the presence of white noise y becomes a stochastic process yt satisfying the 
following stochastic differential equation: 

dy t = y t g(t)dt + edWt (5) 

where g(t) = x(t) — 1, e is a constant with < e fi and W t is the Wiener process. 
Exact solution of (5) is possible [10]: 

yi = G(t,t a )(vo + e^^dW.) (6) 



where 



G(Mo) = e Jt o . (7) 
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The mean value of y t , 

(y t )= G(t,t )y , (8) 

is the solution in the limit e — > 0. The probability distribution of y is Gaussian with 
standard deviation, a y , a function of time given by 

°i = (yt)-(yt) 2 = z 2 G 2 (t, to) J* c^j ds - (9) 

If x = 1 at t = t a then for t - t a > 0(^=) 

al~e 2 ^G 2 (t,t a ) (10) 

where 

= ^7(lW (n) 



2.1. The deterministic limit 

If o y -C (yt) for x > 1 then the noise can be treated as a small perturbation to the 
deterministic solution given by (8) or by 

V = yoe-t™- F <*°» (12) 

where 

= - ln(l + —)- x +- tan" 1 . (13) 

The probability distribution of x max in this case is very narrow and the mean value, 
(s max ), satisfies 

F({x max )) - F(x ) =n(\ny -^ ln(///((x max ))^ . (14) 



2.2. The noise-controlled regime 

If Oy 3> (yt) for x > 1 then it is the noise rather than the initial conditions which 
controls x max . In this case the probability distribution of y for x > 1 is Gaussian 
with negligible mean and exponentially rising standard deviation. Equation (10) 
corresponds to the fact that, for t — t a > (9(-^), the path of any one realisation is 
approximately deterministic (but starting from a level which is random variable). 
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Knowing a y as a function of time, we can calculate the probability that y 2 
is greater than nf(x) at any time. The probability distribution of x max is the 
derivative with respect to x of this probability. The probability that x max lies 
between x and x + dx is therefore R(x)dx where 

q 2 r-Vvf(x) 



R( x ) = 2—^=- / e ^dy 



yJ ~~ (15) 
2 y^fix) -eg? 



27T fif(x) a y 



The maximum value of R(x) is at a: = x where x is defined by the condition 
a/ fif(x) = a y . (Both x and a y are functions of time.) Thus x satisfies 

= e(7ra) 4 (16) 



or equivalent ly 

F(l) - F(x) = n\ lne| + £ ln(///(x)) - j ln(Tra). (17) 

At the very end of the slow phase the simple relationship i = nf{x) breaks down 
because y 2 is no longer negligible. We obtain a more accurate expression for the 
most probable value of x max by replacing x by x c = x — Ax where 

r max 2<it , 1 

Ax = y = -^U. (18) 

We exhibit the probability distribution of x max by defining the random variable 

v as 

Wl^c) 

The probability that t> lies between v and v + cfo is R(v)dv where 

R(v)~-^e v e-^ e2v (20) 

V 27T 



which is the probability distribution of the log of the absolute value of a Gaussian 
random variable with unit variance. Explicit expressions for the mean and variance 
of v exist: 

(<A=_I (7 + l n2 ) (21) 
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where 7 = .577 . . . (Euler's constant) and 



<- 2 >» 2 = y (22) 

so the mean and standard deviation of x max in the noise-controlled regime are given 
by 

<x max )=x c + i^^( 7 + ln2) (23) 

- f/~2 \Z/ 71 Hf(x c ) . , 

Each value of x\ for —1 + 0(y/J2) < x\ < 1 in Figure 2 corresponds to a set of 
initial conditions for the slow phase [1]. The value taken for X2 is the lower of the 
two values given by (14) and (23). When the noise-controlled value (23) is the 
lower, we plot (x max ) + a Xmx and <x max ) - a Xmtx . For the purposes of Figure 2, 
the transition from the noise-controlled to the deterministic regime is sufficiently 
rapid that it is unnecessary to consider the transition region. In the next section, 
however, we derive a more general formula for the probability distribution of x max . 

2.3. A more general formula 

We express the relative magnitudes of the deterministic and noise-driven parts of 
y t via the parameter c defined as 



\Vt\ 



(25) 



which is constant for t — t a > 0(—^) (10). The noise-controlled regime corresponds 
to < c < 1, and the deterministic limit to c > 1. The condition c < 1 for 
t — t a > O(-) can be taken as a test of whether the noise-controlled regime is in 
force. This condition is 

H\ lne| < F(l) - F(x ) + fj, Q ln(Tra) - lnj/ ) • (26) 

A more general formula for R(x) is obtained by allowing the probability dis- 
tribution of y for x > 1 to have non-zero mean. Thus 

(y-(y)) 2 



q 2 rV^Tix) 



dx y^o-y 



R( X ) = —-^I e ^ dy 

-y/JU^) (27) 



' u(e-^ 2 +e-^ 2 ) 9{X) 



2tt V 7 nf(x) 
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Figure 3. The effect of noise on the probability distribution of Xmax- 
Numerically-obtained probability distributions of x ma x with fi = 0.01 and 6 = 1.0. 
The smooth curve in each case is the function R(x) (27) which reduces to (15) in 
the noise-controlled regime. 



where 



u 



(28) 



This probability distribution is compared with numerical results for \x = 0.01 
in Figure 3. In the limit c — > we find the distribution (20). For large c we 
regain the deterministic regime (narrow, Gaussian distribution of x max ). 

Our results are exact for small ft and delta function initial conditions. To 
produce Figure 3 and Figure 4 we take the initial conditions from the corresponding 
deterministic orbit. This gives excellent agreement in the noise-controlled regime, 
where x max is independent of initial conditions, and gives the correct large-c limit 
for (x max ). However the standard deviation cr Xmax is underestimated for nonzero 
c because the probability distribution of x max is carried through the fast phase, 
so the initial conditions for the slow phase vary from cycle to cycle in the noisily 
periodic orbit. The broadening of the probability distribution of x max produced 
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Figure 4. Mean and standard deviation of x m ax as a function of u:\ln e\. 
The curves are predictions obtained from the explicit form for the probability 
distribution of x ma x (27) and some numerical results are shown. The agreement 
between small-u calculations and numerical results is good even for u, as large as 
0.1. The knee in the graph of u Xmax vs u\lne\ corresponds to the end of the 
noise-controlled regime. 

by this is most noticeable in the transition region between the noise-controlled and 
deterministic regimes (e = 10" 16 in Figure 3). 

In Figure 4 we compare numerical results for the most probable value of x max 
and the standard deviation of x max with values calculated using the explicit form 
of the probability distribution (27) . The most probable value of x max we find from 
the (approximate) condition 

u 2 = l + c 2 . (29) 

This corresponds to (%) 2 ) = fif(x) and gives the correct result in the large and small 
c limits. The standard deviation of x max can be written 



(30) 



where f(x) and g(x) are evaluated at the most probable value of x max , h(0) = 1 
and, for large c, h(c) ~ \. The form we have used for h(c) in Figure 4 is /i 2 (c) = 

_ c 2 

Yq^2 + c 2 e~, which we obtained as a fit to numerical integration of (27). 
3. Conclusion 

Slow-fast dynamical systems such as the three-wave resonance system discussed 
here are most conveniently described in terms of a one-dimensional map. Low-level 
white noise has an 0(1) effecttrue which can be calculated by solving a stochastic 
differential equation. In the deterministic limit x max , the turning point of x which 
defines the end of the slow phase, is determined by the initial conditions. In the 
noise-controlled regime, which is in force when /i\ lne| is less than an 0(1) constant 
which depends on initial conditions, x max is a random variable with most probable 
value a function of /x| lne| and standard deviation proportional to fx. 

We know of several other physical contexts which give rise to noise-sensitive 
slow- fast dynamics. One is the shear instability of tall thin convection cells [2, 3]. 
Another is pulsating laser oscillations [11, 12] consisting of short pulses separated 
by long periods of very small intensity. A related problem is that of random pertur- 
bations of heteroclinic attractors [13], where noise controls the length of time spent 
near an unstable fixed point, and a stable homoclinic or heteroclinic orbit provides 
the reinjection. 
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Appendix 

In the three-wave resonance problem considered above, the slow phase is defined by 

x = fif(x), 

(31) 

dy t = g(x)y t dt + edW t 

where f(x) = x + 1- and g(x) = x — 1. We obtain a system often studied in dynamic 
bifurcation theory [14] if we put f(x) = 1 and g(x) = x. Then y t satisfies 

dy t = g(t)y t dt + edW t where g = fit (32) 

Suppose y = yo for some t = to < 0. (In the three-wave resonance problem the 
natural choice is yo = 0(^/Ji) and — /ito = 0(1)0 The next value of g at which 
V = Do? 9i defines a dynamic bifurcation point, occurring later than the 'static' 
bifurcation at g = 0. In fact, in the deterministic limit, g = —g(to). 
The noise-controlled regime is in force if 

tA lne l < ^9 2 (to) + V Q In - In • (33) 
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In this regime, g is a random variable with most probable value, g, given by 

9 2 = 

and standard deviation Og by 



2 = 2//|lne| + ^ln^+2//lny o (34) 

Z 7T 



In this form, our results appear consistent with results obtained numerically [15], 
from an electronic circuit model of a ring laser [16] and with analytical results for 
the laser threshold instability [17]. 
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